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Abstract. There is strong observational evidence that a quasi-thermal 
population of electrons (or pairs) exists in compact X-ray sources. It is, 
^ however, unclear what mechanism thermalizes the particles. Here, two 

qj , processes. Coulomb scattering and synchrotron self-absorption, that may 

pLn ' be responsible for the thermalization, are reviewed. The parameter spaces 

in which respective process dominates are given. 

While the Coulomb thermalization mechanism is well-known, this is 
not the case for the synchrotron self-absorption thermalization. We give 
the arguments that synchrotron self-absorption must act as a thermal- 
izing mechanism in sufficiently compact sources. The emitting and ab- 
(^ I sorbing electrons then exchange energy efficiently with the self-absorbed 

CN ■ synchrotron radiation field and are driven towards a relativistic or mildly 

CN I relativistic thermal distribution in a few synchrotron cooling times (the 

2r ' "synchrotron boiler"). 

o 

-4—' ' Observations with hard X-ray /7-ray satellites such as CGRO OSSE, RXTE, and 

j^ I BeppoSAX indicate that the X/7 spectra cut off at a few hundred keV for the 

majority of active galactic nuclei (see Zdziarski, this volume; Matt, this volume) 

and for the hard states of galactic black hole candidates (see Zdziarski, this 

^ ' volume; Grove, this volume). It is generally believed that Comptonization by a 

quasi-thermal population of electrons (or pairs) is responsible for the formation 
of the X/7 spectra from these sources. In the spectral modeling codes (e.g., 
Poutanen & Svensson 1996), it is normally assumed that the Comptonizing par- 
ticles have a Maxwellian distribution. It has been pointed out several times that 
thermalization by Coulomb scattering may not be fast enough as compared to 
various cooling mechanisms (such as Compton cooling) and that the particle 
distribution therefore will differ from a Maxwellian (e.g., Dermer &: Liang 1989; 
Fabian 1994). On the other hand, it has been noticed that another thermaliza- 
tion mechanism, synchrotron self-absorption, may operate in compact plasmas 
(Ghisellini, Guilbert, & Svensson 1988; Ghisellini &: Svensson 1989). 

Here, we review the physics of these two thermalization mechanisms, and 
explore in which contexts each of them may operate. 
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Introduction 



2. Thermalization by Coulomb Scattering 

The approximate time scale, tc, for thermalization by Coulomb (M0ller) scat- 
tering between electrons in nonrelativistic plasmas has been known since long 
(e.g., Spitzer 1956; see Stepney 1983 for relativistic corrections): 

tc = 4 JZ-e3/2(vri/2 _ 1.291/4 + 26^2), (i) 

m A 

where tx = (^e'^Tc)^^ is the Thomson time, ng is the electron density, InA ~ 
10 — 20 is the Coulomb logarithm, and & = kTe/nieC^ is the dimensionless tem- 
perature. Only with the advent of X-ray astronomy emerged the understanding 
that the electrons may reach mildly relativistic temperatures of 10^ K or more, 
and that electron-positron pairs may be created. In such conditions, electron- 
positron (Bhabha) scattering also becomes important. Here, the term "Coulomb 
scattering" will be used for all types of "Coulomb" interactions. 

2.1. The Relaxation Process 

Detailed numerical simulations of the relaxation and thermalization process were 
performed by Dermer & Liang (1989). As small angle scatterings normally 
dominate and the fractional energy change per scattering is small, a Fokker- 
Planck approach is appropriate. Dermer & Liang evaluated the energy exchange 
and diffusion coefficients assuming the plasma to have a Maxwellian distribution. 
Using these coefficients in the simulations is equivalent to studying the relaxation 
of a test particle distribution in a Maxwellian background plasma. Nayakshin &: 
Melia (1998) relaxed the Maxwellian assumption and computed self-consistent 
Fokker-Planck coefficients using the real particle distributions. A Monte Carlo 
approach was taken by Pilla &: Shaham (1997) who treated the time evolution 
of both the pair and photon distributions in an infinite system and who besides 
Coulomb interactions, bremsstrahlung, and Compton scatterings also included 
pair production and annihilation. 

Figure 1 (from Dermer & Liang 1989) shows the relaxation of a test electron 
distribution in a background thermal electron plasma of temperature 511 keV, 
or, equivalently, = 1. The test electrons had initially a Gaussian distribution 
centered at 1 MeV and a FWHM of 0.28 MeV. The diffusion process dominates 
initially and broadens the electron distribution that first relaxes at lower ener- 
gies and only later the Maxwellian tail forms. Figure 7 in Nayakshin & Melia 
(1998) shows the similar process but here all electrons are initially Gaussian 
(i.e., there is no Maxwellian background). If the initial distribution is broader 
than a Maxwellian then the energy exchange coefficient dominates initially and 
the low energy electrons gain energy while the higher energy electrons loose 
energy thereby narrowing the distribution towards a Maxwellian (see fig. 6 in 
Nayakshin & Meha 1998). 

2.2. Influence of Cooling Processes on the Steady Electron Distri- 
bution 

With increasing temperature, the Coulomb energy exchange rate decreases. 
Various cooling processes such as bremsstrahlung, Compton cooling and syn- 
chrotron cooling increases with temperature and eventually the thermalization 
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Figure 1. Thermalization through electron-electron scattering of an 
initially Gaussian electron test distribution (curve a) centered at 1 
MeV and a FWHM of 0.28 MeV in a background thermal electron 
plasma having a Maxwellian distribution of temperature 511 keV, i.e., 
= 1 (curve b). The dashed curves show the relaxing distribution 
at different times, assuming an electron background density of 1 cm~'^. 
The thermalization time scale from eq. (1) is 2.6 x 10^^ s. From Dermer 
k Liang (1989). 



process will be inhibited by the cooling first noticeable as a truncation of the 
Maxwellian tail. Stepney (1983) noticed that bremsstrahlung cooling will pre- 
vent thermalization for temperatures larger than about 5 xlO^° K, and Baring 
(1987) performed further analysis for additional cooling processes, as did Ghis- 
ellini, Haardt, & Fabian (1993). 

Including Compton and synchrotron losses in the Fokker-Planck equation 
allows for the determination of the steady distribution function under the influ- 
ence of these cooling processes. Results from Dermer & Liang (1989) are shown 
in Figure 2. It is seen that, for increasing energy densities of radiation and mag- 
netic fields, the high energy tail of the electron distribution becomes increasingly 
truncated and the effective temperature of the distribution becomes smaller. 

What are then the conditions for losses to dominate over Coulomb thermal- 
ization (see, e.g., Fabian 1994)? The nonrelativistic cooling time scale can be 
written as (see Coppi, this volume) 

tcool-R/[ciB{l + Ur,d/UB)], (2) 

where R is the size of the region, Ib is the magnetic compactness defined in 
equation (4) below, and C/rad, Ub are the energy densities in radiation and mag- 
netic fields, respectively. Comparing with equation (1), one finds that Coulomb 
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Figure 2. Truncation of a thermal electron distribution with tem- 
perature 511 keV due to Compton and synchrotron losses. Solid curve 
shows a Maxwellian of temperature 511 keV. Dashed curves show the 
relaxed distributions at different values of the parameter (?7rad + t^B)/^e 
in units of MeV. From Dermer & Liang (1989). 



scattering cannot maintain a Maxwellian when 



^B(l + ^rad/f/B)> 



where tt = ngCT-R is the Thomson depth. 
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2.3. Can Coulomb Collisions Thermalize Pair Coronae and Active 
Pair Regions? 

Two different numerical codes (Stern et al. 1995a; Poutanen & Svensson 1996) 
have been used to study radiative transfer and Comptonization in pure pair 
coronae in energy and pair balance (see, e.g., Svensson 1997 for a review). For 
coronae of a given geometry and in energy balance, there exists a unique Tg — 
Tt relation, where Tg is the volume-averaged coronal temperature and tt is a 
characteristic Thomson scattering optical depth of the coronal region. In Figure 
3a, this relation is shown for different geometries. The results for active regions 
are connected by dotted curves. For comparison we also show the slab results 
from Stern et al. (1995b) using an iterative scattering method code {dashed 
curve) . 

Solving the pair balance for the obtained combinations of (G, tt) gives a 
unique dissipation compactness, £diss (see Ghisellini & Haardt 1994 for a dis- 
cussion). Here, the local dissipation compactness, ^diss = {Ldiss/h){aT/meC^), 



characterizes the dissipation with -L^iss being the power providing uniform heat- 
ing in a cubic volume of size h in the case of a slab of height h, or in the whole 
volume in the case of an active region of size h. Figure 3b shows the Q vs. ^diss 
relations. 
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Figure 3. Dimensionless volume-averaged temperature, & = kTe/nieC^, 
vs. Thomson scattering optical depth, tt, in panel (a), and vs. dissipation 
compactness, ^diss = (Ldiss/h) {aT/nieC^) in panel (b), for a steady X-ray 
emitting plasma region in pair and energy balance on or above a cold disk 
surface of black body temperature, fcTbb = 5 eV. The plasma Compton scat- 
ters reprocessed soft black body photons from the cold disk surface. Solid 
rectangles and dashed curve show results from a nonlinear Monte Carlo code 
(Stern et al. 1995a) and an iterative scattering method code (Poutanen, & 
Svensson 1996), respectively, for the case of a plane-parallel slab corona. Re- 
sults using the Monte Carlo code for individual active pair regions are shown 
for hemispheres located on the disk surface; surface spheres also located on the 
surface {underlined spheres); spheres located at a height oi0.5h (spheres), 1ft. 
(diamond), and 2h (triangle), where h is the radius of the sphere. The results 
for each type of active region are connected by dotted curves. The dash-dotted 
and dash- dot- dot- dotted curves in panel (b) show the critical compactness as 
function of O above which thermalization by M0ller and Bhabha scattering 
is not achieved for the cases of pair slabs and surface spheres, respectively. 
(Unpublished results by Stern et al.; see also Stern et al. 1995b; Svensson 
1997). 



The question arises whether the electrons can thermalize or not for the 
conditions, 0, tt, and ^dissj in Figure 3. Energy exchange and thermalization 
through M0ller (e^e^) and Bhabha (e^e^) scattering compete with various loss 
mechanisms, with Compton losses being the most important for our conditions. 
The thermalization is slowest and the Compton losses largest for the higher en- 
ergy particles in the Maxwellian tail. Instead of using the approximate equation 
(3), we use the detailed simulations by Dermer & Liang (1989, their fig. 8) to find 



the critical compactness above which the deviation of the electron distribution 
at the Maxwellian mean energy is more than a factor e ~ 2.7. The dash-dotted 
and dash- dot- dot- dotted curves in Figure ^ show this critical compactness for 
slabs and for surface spheres, respectively. In agreement with Ghisellini, Haardt, 
& Fabian (1993), we find that M0ller and Bhabha scattering cannot compete 
with Compton losses in our pair slabs and active regions. 

The problem then arises of what mechanism can thermalize the apparently 
thermal electron distribution in compact sources. One such mechanism is cy- 
clo/synchrotron absorption. 

3. Thermalization by Cyclo/Synchrotron Absorption 

3.1. A Brief History of Synchrotron Thermalization 

Ever since the classical interpretation by Shklovskii in the 1950s of the radiation 
from the Crab nebula as being synchrotron radiation, this process has played an 
important role in our interpretation of the non-thermal radiation from a wide 
variety of astronomical objects. In general, the electron distribution has been 
assumed to be a power law or nearly a power law. Much less attention has 
been paid to what happens to the electron distribution at self-absorbed electron 
energies. 

The theory for synchrotron radiation was developed in the 1950s (see, e.g., 
reviews by Ginzburg and Syrovatskii 1965, 1969; Pacholczyk 1970). The emis- 
sion and absorption coefficients for single relativistic electrons as well as for en- 
sembles of relativistic electrons having power law or thermal distributions were 
calculated in the 1950s and 1960s. For power law distributions, the absorption 
coefficient increases towards lower photon energies. Below some photon energy, 
fabs! the source becomes optically thick to synchrotron self- absorption resulting 
in an intensity proportional to i/^'"^, obtained from the ratio of emission and ab- 
sorption coefficients (Le Roux 1961). A finite i^abs> of course, requires that the 
source is finite. Below we call this and its consequences "finite source effects". 

This early work was, in general, applied to extended sources where the 
cooling time at self-absorbed particle energies is longer than other relevant time 
scales (such as the age of the source or the dynamical time scales). It was 
therefore natural to assume that the self-absorbing electron distribution below 
the Lorentz factor, 7abs) of the electron (emitting at the frequency t'abs where 
the source becomes optically thick) is unaffected by self-absorption and simply 
maintains the power law distribution of the injected electrons. 

In the late 1960s, it became increasingly clear that, in compact sources or on 
long time scales, the self-absorbed electron distribution A^(7) will evolve under 
the influence of synchrotron emission and absorption. What are then the possi- 
ble equilibrium solutions at self-absorbing Lorentz factors towards which A^(7) 
would relax? In the important papers by Rees (1967) and McCray(1969), it was 
shown that power law distributions, A'^(7) c>c 7"^^ with s = 2 and 3, are equilib- 
rium solutions to the kinetic equations. Rees (1967), however, also found that 
the solution with s = 3 is unstable and would evolve away from s = 3 if slightly 
perturbed. McCray (1969) showed this explicitly by numerically calculating 
the time dependent evolution of initial power law distributions in an infinite 
source. Rees predicted and McCray confirmed that the high energy electrons 
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in a flat (s < 3) initial power law would tend to evolve into a quasi-Maxwellian 
distribution. McCray (1969), furthermore, emphasized the importance of finite 
source effects on the evolution. In particular, for power laws with s < 3, the 
self-absorbing electrons would gain energy absorbing slightly more energy than 
they emit, while the electrons radiating in the optically thin limit lose energy 
by radiating much more energy than they absorb. All electrons would therefore 
tend to gather at 7abs developing a peak there (as was already emphasized by 
Rees 1967). 

It must be emphasized that the relaxation of self-absorbing electrons takes 
place through the energy exchange with the radiation field, which in its turn 
is determined by the particle distribution. This is the "synchrotron boiler", a 
terminology coined by Ghisellini, Guilbert, &: Svensson (1988). 

3.2. Rise and Fall of yet another Paradigm 

In a series of papers in the 1970s, Norman and coworkers further developed the 
concept Plasma Turbulent Reactor (PTR) introduced by Kaplan and Tsytovich 
(1973). Originally the turbulence feeding the electrons was thought to be plas- 
mons. In practice, however, the PTR is exactly the self-absorbing synchrotron 
source considered here, as photons are the only plasma modes with sufficiently 
small damping rate that they can mediate energy transfer from one electron 
to another (Norman 1977). Norman and ter Haar (1975) and Norman (1977) 
essentially repeated the analysis of McCray (1969) using quite different nota- 
tion and definitions but arriving at the same conclusions that A^(7) oc 7^^ and 
A^(7) oc 7~^ are the only steady power law equilibrium solutions. It is important 
that they noted that the A^(7) oc 7"^ solution corresponds to a finite electron 
flux upwards along the energy axis, while A^(7) oc 7"^ corresponds to zero elec- 
tron flux. They argued that A^(7) oc 7^^ was the most physical solution as the 
synchrotron time scales establishing this distribution are shorter than other time 
scales. Although being aware of possible finite source effects, they considered 
them not to influence the electron distribution at Lorentz factors <ti 7abs- 

The self-absorbed solution, N{'y) oc 7~^, was considered sufficiently impor- 
tant in explaining power law spectra from a variety of sources that Norman 
and ter Haar (1975) called the PTR a new astrophysical paradigm. Norman 
and coworkers, however, do not seem to have considered the stability of the 
A^(7) oc 7~^ solution. 

The work of Rees (1967) and McCray (1969) indicates that a Maxwellian 
distribution may be the only stable equilibrium solution. This was, however, not 
rigorously established causing Ghisellini, Guilbert & Svensson (1988, GGS88) 
to numerically determine the steady solutions of the kinetic equations including 
physical boundary conditions (i.e., correct Fokker-Planck coefficients at subrel- 
ativistic energies, and the accounting for finite source effects at large energies). 
As 7abs typically is of the order 10-100 in compact radio sources and the develop- 
ment of the self-absorbed distribution takes place at mildly relativistic energies, 
they used expressions and equations valid at any energy. Furthermore, in or- 
der to obtain steady solutions the particle injection had to be balanced by a 
sink term (escape or reacceleration) . Injecting a power law proportional to 7"'^ 
(i.e., with the equilibrium slope) GGS88 found that the steady solution was a 
Maxwellian with a temperature corresponding to the mean energy of the in- 



jected electrons. Similarly, an injected power law proportional to 7^ (essentially 
corresponding to monoenergetic injection at some large Lorentz factor ^ 7abs) 
led to the establishment of a Maxwellian just below 7abs- The injected electrons 
cool until reaching 7abs where they thermalize exchanging energy with the self- 
absorbed radiation field. Additional studies were made by de Kool, Begelman, 
& Sikora (1989) and Coppi (1990). With these works it appears that the PTR 
paradigm of Norman and ter Haar (1975) has been shown to be invalid. 

3.3. Relaxation by Cyclo/Synchrotron Absorption 

The works so far that explicitly have demonstrated the formation of a Maxwellian 
through synchrotron self-absorption are the numerical simulations of GGS88, 
Coppi (1990), and Ghisellini, Haardt, & Svensson (1998). Here, we review some 
of the results in the last paper. In the numerical simulations, a kinetic equation 
for the electron distribution is solved. The kinetic equation, which also in this 
case takes the form of a Fokker-Planck equation (see derivation in McCray 1969), 
includes Compton and synchrotron cooling, synchrotron absorption (heating), 
electron injection, and electron escape. Even though various source geometries 
are discussed, the radiation field is assumed to be given by the steady slab 
solution, which is correct to the order of unity. 

The simulations consider a region of size R with a magnetic field of strength 
B, into which some distribution of electrons are injected with a power L. In a 
steady state, this power emerges either as radiation or as the power of escaping 
electrons. The electrons are assumed to escape at the speed fesc = c/3esc = R/icsc 
where tesc is the escape time. Convenient parameters describing compact sources 
are the injection compactness, £inj, and the magnetic compactness, 1-q-, defined 
as 

^-j = ^ ^ ; Ib = ^RUb. (4) 

where ut is the Thomson cross section and Ub = B'^/Stt is the magnetic field 
energy density. Note that C/rad/t^B ~ (9/167r)(£inj/^B)(l + 7"T), where the numer- 
ical factor is dependent on the source geometry. The problem we consider has 
the following parameters: ^inj, ^B; /3esc; and either R or B. Further parameters 
are those describing the shape of the injected electron distribution. 

For steady state electrons emitting and absorbing synchrotron photons, the 
cooling (emission) and absorption/diffusion time scales are balanced and thus 
equal. The synchrotron cooling time scale can thus be taken to be the ther- 
malization time scale. From Equation (2), it is then clear that self-absorbing 
electrons will thermalize before they escape when £b ~ 1- The simulations shown 
below have £b = 10 and 30. 

Figure 4 shows the relaxation due to cyclo/synchrotron absorption of an 
electron distribution towards the equilibrium Maxwellian distribution. The in- 
jected electrons have a Gaussian energy distribution peaking at 7 = 10. Each 
curve is labeled by the time in units of R/c. The shape of the equilibrium dis- 
tribution is reached in about ^ O.lR/c, about equal to the cyclo/synchrotron 
cooling time. With the assumed input parameters, the synchrotron terms (emis- 
sion, absorption and energy diffusion) in the kinetic equation are dominant over 
Compton losses. Gains and losses in this case almost perfectly balance. As a 
result the equilibrium electron distribution is a Maxwellian. Figure 4 also shows 
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Figure 4. Electron distribution, r(p) = a^RN^p), evolving due to 
cyclo/synchrotron eniission/absorption/diffusion. Curves are labeled 
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and /3esc = 1- The injected distribution is 
= 10. From Ghisellini, Haardt, & Svensson 



that the high energy part of the Maxwellian distribution is formed earlier than 
the low energy part, due to the higher efficiency of exchanging photons of the 
high energy electrons. A slower evolution takes place after 0.1 (i?/c), as the bal- 
ance between electron injection and electron escape is achieved on a time scale 
of a few tesc- Only then have both the shape and the amplitude of the electron 
distribution reached their equilibrium values. 

3.4. Influence of Cooling Processes on the Steady Electron Distri- 
bution 

The equilibrium distributions for different values of the injected compactness are 
shown in Figure 5. The magnetic compactness is set to Ib = 30, corresponding to 
B = W^ G for R = W^^ cm (from Eq. 4). In all cases, the injected distribution 
is a peaked function with an exponential high energy cut-off. The mean injected 
Lorentz factor is < 7 >~ 5 and essentially all electrons are below 7abs- It is ap- 
parent from Figure 5 that the electron distribution is a quasi-Maxwellian at all 
energies as long as £inj <C Ib ■ This is a consequence of an almost perfect balance 
between synchrotron gains (absorption) and losses, while Compton losses are 
only a small perturbation. As ^mj increases towards Ib, Compton losses become 
increasingly important, competing with the synchrotron processes. At high en- 
ergies, losses overcome gains, and the electrons diffuse downwards in energy, 
until subrelativistic energies are reached. In this energy regime, the increased 
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Figure 5. Steady equilibrium electron distributions due to cy- 
clo/synchrotron eniission/absorption/diffusion and Compton cooling 
for different injected compactnesses (decreasing from top to bottom). 
Further parameters are ^b = 30, R = 10^^ cm, /3esc = 1- The corre- 
sponding magnetic field is i3 = 10^ G. The mean injected Lorenz factor 
is about 5. Increasing £inj implies increased Compton cooling resulting 
in a shift of the quasi-Maxwellian towards lower temperatures. From 
Ghisellini, Haardt, & Svensson (1998). 



efficiency of synchrotron gains (relative to losses) halts the systematic downward 
diffusion in energy, and a Maxwellian can form (see Ghisellini & Svensson 1989). 
The temperature of this part of iV(7) can be obtained by fitting a Maxwellian 
to the the low energy part of the distribution, up to energies just above the peak 
of the electron distribution. 

The resulting temperatures are plotted in Figure 6 as a function of ^inj. 
For £inj ;$ 1, the temperature is approximately constant, while it decreases for 
£inj ~ 1. From Equation (3) (with ?7rad/f^B set to zero), we see that thermal- 
ization by synchrotron self-absorption dominates when £3 > "^T In A/40^/^ (as- 
suming 0^1). The Coulomb process thus dominates for small temperatures 
and large tt (i.e., large electron densities). We need to know tt for our simu- 
lations. The balance between electron injection and escape in our model gives 
a Thomson optical depth of tt = (3/47r) (£inj //3esc < 7 >)• For the simulations 
in Figure 5, the optical depth increases from tt = 5 x 10~^ for t\m = 0.1 to 
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5 for 



'inj 



100. Using the expression for tt, we find that thermal- 



ization by synchrotron self-absorption then dominates over Coulomb scattering 
for temperatures > 0.11(lnA/ < 7 >)^/'^(Anj/^B)^''^- For the parameters 
of the simulations in Figures 5 and 6 and InA = 20, the condition becomes 
> 0.03(^inj)^'^, which is plotted as the solid line in Figure 3. One sees that 
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Figure 6. Temperatures of the quasi-Maxwellian part of steady equi- 
librium electron distributions shown in Fig. 5. Above the solid line, 
synchrotron self-absorption dominates over Coulomb exchange as the 
thermalization mechanism. The solid dots show the Comptonization 
y-parameter. From Ghisellini, Haardt, & Svensson (1998). 



that synchrotron self-absorption dominates the thermalization for all cases with 
^inj smaller than about 10 . For the cases £inj = 30 and 100, one cannot neglect 
Coulomb thermalization. 



3.5. Spectra from Steady Electron Distributions 

In Figure 7, the radiation spectra corresponding to four of the equilibrium elec- 
tron distributions in Figure 5 are shown. Each spectrum consists of several 
continuum components: 

• a self-absorbed synchrotron spectrum (S); 

• a Comptonized synchrotron spectrum (SSC); 

• a reprocessed thermal soft component (bump); 

• a component from Comptonization of thermal bump photons (IC); 

• a Compton reflection component. 

Details of the spectral calculations are given in Ghisellini, Haardt, & Svens- 
son (1998). Some features in Figure 7 may be noticed. For £inj < 1, the Compton 
y-parameter is less than unity (see Fig. 6) making the Compton losses relatively 
unimportant relative the self-absorbed synchrotron radiation. The large value of 
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Figure 7. Radiation spectra calculated using four of the electron dis- 
tributions shown in Fig. 5. The different spectral components are 
the reprocessed thermal bump and IC components (dashed curves), 
the synchrotron and SSC components (dash-dotted curves), and the 
Compton reflection component (dotted curves). The total spectra are 
shown by the solid curves. Note that the chosen parameters {R = 10^^ 
cm, and the corresponding magnetic field B = 10^ G) correspond to 
conditions expected near black holes in active galactic nuclei. From 
Ghisellini, Haardt, & Svensson (1998). 



B makes the Comptonized spectra bumpy. The 2-10 keV band is dominated by 
the SSC component, rather than by the IC. The thermal bump and the X~ray 
flux are thus not directly related. This is contrary to the common interpreta- 
tion of the X-ray emission in Seyfert galaxies as being due to Comptonization 
of thermal bump photons. 

For ^inj > 1, the Compton cooling dominates and limits the y-parameter to 
unity. The smooth IC power law dominates over the S and SSC components. 
For ^inj ^ 3, the high energy spectral cut-off can be described by an exponential, 
since the electron distribution is a quasi-Maxwellian in the entire energy range. 
For £inj ^ 3, the electron distribution is more complex (see Fig. 5), resulting in 
a more complex spectral cut-off. 

The choice of i? = 10^^ cm (or B = 10^ G) in Figure 7 corresponds to the 
case of active galactic nuclei (AGN). Ghisellini, Haardt, & Svensson (1998) also 
study the case of galactic black holes choosing R = 10^ cm (or B = 10^ G) for 
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the same set of i^ and ^inj. The large magnetic field and small size move the 
S and bump peaks to larger frequencies. The main difference is that now the 
synchrotron component is not completely self-absorbed leading to an optically 
thin synchrotron component from the highest energy electrons. More soft syn- 
chrotron photons enhances the SSC component relative the IC component as 
compared to the AGN case. 

4. Final Remarks 

First, we note that £b > 1 is needed for synchrotron self-absorption to operate 
efficiently (from eq. 2). A Maxwellian is then formed with the same mean energy 
as the injected electrons, assuming that self-absorption operates at essentially 
all electron energies of interest. If, furthermore, Compton cooling is important, 
i.e., if Anj > ^B) then the Maxwellian is modified and is shifted to lower energies 
(temperature). 

Second, we note that the criterion for Coulomb vs synchrotron thermaliza- 
tion in the case when all electrons are self- absorbing (i.e. radiating optically 
thick synchrotron radiation) is more complex (eq. 3). Essentially, the same cri- 
terion is valid for comparing Coulomb thermalization in the case when the major 
part of the electrons radiate optically thin synchrotron radiation (or Compton 
radiation) truncating the Maxwellian. Which of the two cases that is valid de- 
pends on whether the energy, (7abs — ^)meC^, of the electrons radiating at the 
photon energy where the absorption optical depth is unity is much larger or 
much smaller than fcTg. 

There should also be a region in parameter space where both Coulomb and 
synchrotron thermalization operates simultaneously (assuming that most elec- 
trons radiate optically thick radiation). Here, Coulomb thermalization should 
dominate at lower electron energies and synchrotron thermalization at larger 
energies. However, nobody seems so far to have solved the Fokker-Planck equa- 
tions to study the thermalization process including both Coulomb scattering 
and synchrotron self-absorption. 

Ultimately, the thermalization process should be put in a realistic con- 
text. Mahadevan & Quataert (1997) studied the importance of thermalization 
in advection-dominated flows onto black holes under the conditions considered 
in such flows (e.g., close to free-fall, equipartition magnetic fields). Comparing 
the thermalization time scales with the accretion time scale (equivalent to tcsc in 
our discussion), they found that thermalization did not occur at large radii and 
small accretion rates. However, at sufficiently large accretion rate, synchrotron 
thermalization becomes important, and at even larger rates (and thus larger 
densities) Coulomb thermalization starts operating. 

Another scenario for generating the X-ray radiation from compact sources 
is that of a corona or magnetic flares atop an accretion disks. The typical condi- 
tion for the flare regions is that the magnetic energy density should dominate the 
radiation energy density, i.e., that Ib ~ ^inj > 1, which should ensure that cy- 
clo/synchrotron self-absorption acts as a very efficient thermalizing mechanism 
in such regions. 
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